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THE  CHORD  DISTRIBUTION  FOR  A 
RIGHT  CIRCULAR  CYLINDER 
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In  the  evaluation  of  the  effects  of  an  omnidirectional  fluence  on  a  sensitive  volume  the  chord 
distribution  is  an  appropriate  tool.  It  expresses  the  frequency  of  occurrence  of  chords  of  a  given 
length  for  a  known  convex  volume  (in  the  three  dimensional  case).  In  particular  in  single  event  upset 
calculations  one  often  calculates  the  upset  rate  of  an  omnidirectional  fluence  passing  through  the  sen¬ 
sitive  volume  of  a  memory  cell  by  using  the  cosmic  ray  approximation,  which  makes  use  of  the  chord 
distribution  of  the  sensitive  volume. 

This  report,  after  reducing  the  differential  chord  distribution  for  a  right  circular  cylinder  to 
Carlson’s  elliptic  integrals  and  standard  integrals,  exhibits  the  evaluation  program  and  numerical  and 
graphical  results.  Various  checks  support  the  validity  of  the  evaluation  and  precision  is  discussed. 
The  chord  distribution  is  found  to  exhibit  a  cusp  and  a  jump.  A  plausibility  argument  is  offered  for 
the  existence  of  these  features  and  an  analysis  of  consequences  in  signal  event  upset  calculations  is 
made. 

In  Reference  1,  Kellerer  points  out  some  basic  distinctions  in  probability  measure  of  uniform 
distributions  in  two  and  three  dimensions,  establishes  some  background  and  derives:  1)  the  path 
length  distribution  in  a  convex  body  for  finite  track  lengths  and  2)  the  chord  distribution  in  general 
cylinders  given  the  chord  distribution  for  the  cross-section.  The  present  interest  is  only  in  the  circular 
case  of  tne  latter  and  further  in  the  density  distribution,  Kellerer' s  formula  56: 


c(s)  = 


where 


7 x{d  +  2  h ) 


h  “  oo 

-s  xif(sx)  +  2x2F(sx)  dx  H — -  H(s-h)\  F(x)dx  [  (1) 


F(t)  =  ■vT-  t2/d2 ,  fit)  =  t/d'Td2  -  t2,a  = 


H(s  -  h)  vi  -  h2/s 2 


and  H  is  the  step  function  iH(x)  =  1  for  x  >  0,  H(x)  =  0  for  x  <  0).  Here  d  is  the  diameter  of 
the  base,  h  is  the  height,  and  s  is  the  chord  length.  Also  this  expression  is  normalized.  Kellerer  has 
given  here  the  expression  for  a  general  convex  base  which  has  a  known  chord  distribution  F  (or 
known  /,  the  latter  being  the  density  distribution  and  /  =  —F').  Note  F,  f  =  0  for  t  >  d.  (In 
the  theory  of  generalized  functions  one  regularizes  these  conditions  at  t  =  d  by  writing  F  —  HF  so 
/  —  Hf  +  6F ,  5  being  the  delta  function  so  one  has  H'(x )  =  5(.v).)  Writing  one  expression  and 
separating  into  terms  indexed  by  the  power  of  h  ,  one  has 


c(s)  = 


ird{d  +  2  h) 


!(l  -  x2)(d2  -  s2x2) 


~  c .  +  cn  +  r, 
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The  last  two  terms  are  standard  integrals  but  the  first  should  be  recognized  as  an  elliptic  integral. 
One  notes  the  maximum  chord  exceeds  d  and  that  at  s  =  d  at  the  upper  limit  the  two  square  root 
poles  coalesce  to  form  what  would  appear  to  be  a  non-integrable  singularity.  That  this  is  only  appear¬ 
ance  may  be  seen  by  applying  the  delta  function  substitution  in  the  above  parenthetical  note.  Even  if 
one  were  inclined  to  try  direct  numerical  evaluation,  he  would  not  be  surprised  when  numerical  prob¬ 
lems  arose  while  trying  to  evaluate  near  this  pole.  The  fact  that  this  feature  is  a  standard  one  for  an 
elliptic  integral  should  be  sufficient  incentive  to  pay  the  penalty  of  the  algebra  necessary  to  reduce  this 
integral  to  standard  elliptic  integrals.  After  this  penalty  is  paid  the  pole  causes  no  further  problem. 

Still  this  algebra  is  sufficiently  onerous  that  it  provides  the  main  incentive  for  this  report.  Hav¬ 
ing  once  been  forced  to  do  this  work,  one  has  good  reason  not  to  lose  the  technique.  In  addition  we 
were  forced  by  the  circumstance  of  having  only  a  single  software  source  [2]  to  use  the  modem  listing 
of  standard  elliptic  integrals  advanced  by  Carlson  [3]  rather  than  the  Legendre  standard  forms.  This 
was  fortunate  as  Carlson’s  method  takes  full  advantage  of  the  considerable  symmetries  in  elliptic 
integrals,  resulting  in  certain  evaluation  advantages. 

In  cb  s  =  d  is  a  branch  point  requiring  a  change  in  form.  This  is  done  by  the  change  of 
integration  variable,  y  =  sx Id.  The  resulting  integral  has  the  same  form  provided  one  changes  the 
parameter’s  definition  at  the  branch  point:  m  =  s2/d2  for  s  <  d ,  m  =  d2/s2  for  s  >  d .  Then  for 
5  >  d  one  obtains 


Ci(j)  = 


ird(d  +  2 h)  }*a/d  V(i  _  my2)(i  -  y2) 


Thus  both  branches  of  this  elliptic  integral  can  be  reduced  using  the  same  integral  form.  To  aid  the 

.1  .1  .a 

reduction,  use  of  a  zero  limit  may  be  achieved  by  j  ~  J0  _  J0-  In  following  the  integration 
variable  is  transformed  by  y  =  -x2, 


n  t  ,  r“  x4dx  1  f° 


(~yf/2dy 
(1  +>’)(!  +my ) 


Then  using  SFAM  formula  8.1-1  with  the  substitutions  a  =  1,  a'  =  5/2,  y  =  0,  x  =  - u 2,  z,  =  1. 
i  =  1,2,3,  W|  =  1 ,  w2  =  m ,  w3  —  0,  b  [  2  =  1/2,  i>3  =  5/2,  one  has 


Rr(u,m)  =  j  B  1,  u5R  5  y,  y,  y,  1  -  u2,  1  -  mu2,  1 

L  J  ;  l 


2 


since  the  beta  function  fl,  evaluates  to  2/5.  The  R  function  is  related  to  the  hypergeometric  series 
and  is  used  by  Carlson  to  present  a  unified  treatment  of  many  special  functions.  This  function  is 
homogeneous:  simultaneous  permutation  of  the  first  three  and  last  three  arguments  leaves  it 
unchanged.  The  local  software  librar>  [2]  includes  the  following  standard  elliptic  integrals 


Rf(z)  =  ~  JQ  dt  [(r  +  z,)(r  +  z2)(r  +  z3)J 


-t  2  ’  2  ’  2 


*D<z>  =  |  Ot[(t  +  Zl)('  +  Z2)U  +  Zi)i]  2  =  2'  2'2'z 
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where  z  is  understood  as  a  vector,  z  =  (Z\,  z2,  Z3)  and  the  integral  representations  are  from  SFAM 
formula  8.3-1,  The  intention  then  is  to  reduce  Ry  to  RF ,  RD  and  polynomials  following  SFAM  para¬ 
graph  9.3.  Applying  the  formula  of  SFAM  exercise  5.9-8  with  the  substitutions 
t  =  -3/2,  c  =  5/2,  h'i  =  3/5,  one  has 


-(Zi  -  z2)(zi  ~  z3)  R  5 


All 
2'  2’  2  ’  Z 


~  Z2Z3  R  _5. 
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1  1 


2  ’  2  ’  2'Z 


+ 

r  > 

2 

Z\  -  —  (Z2  +  z 3) 
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The  second  term  is  an  RD  but  the  other  terms  need  further  reduction.  For  the  first  term  one  notes 
c  =  a  so  a  ’  =  0  and  applies  SFAM  formula  6.6-5  to  obtain 


R-± 

2 


3  1  1 


2  ’  2  ’  2  ’  Z 


=  (zi3z2 


(9) 


For  the  third  term  apply  SFAM  formula  5.9-7  with  t  =  -1/2,  b,  =  1/2.  c  =  3/2  and  solve  for  the 


R  1  term  to  obtain 


R 


3  1  1 


2'  2  ’  2  ’Z 


=  j  Rf(z)  ~  ys,/?D(z) 


(10) 


Then  substituting  these  results  into  (8)  with  z  1  =  1 ,  z2  =  1  -  m2,  z3  =  1  -  mu2,  and  then  into  (5), 
one  has 


Ry(u,m)  =  — — 
im 


Ml  -  «2)(1  -  mu2)  +  —  (m  +  1  )u2RD(z)  ~  RF(z) 


(ID 


When  substituting  from  (6)  and  (7)  the  homogeneity  property  has  been  used  to  permute  arguments. 


In  case  u  =  1  here  the  elliptic  integrals  are  complete.  Defining  Rx  by  — —  Rx  =  Ry,  one  has  for 

3  m 


s  <  d 


c,(.v)  = 


8/1 


ir(<7  +  2 h  )s  3 


~  Ry  (m  ,  m ) 


u  - 1 


(12) 


and  for  s  ~>  d 


C,(J)  = 


«/!(/ 


1  ■/  1 


7 r(d  +  2h  )s2  3 
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—  Rx  <u.  m) 


(13) 


U  —  S(2  Id 


One  notes  that  if  5  <  h.  a  =0  and  the  lower  limit  term  vanishes  in  both.  This  condition  is  a 
second  branch  point  in  the  range  of  .v . 
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The  last  term  is  integrable  by  standard  methods 


c20s)  = 


4  h2d 


H(s  -  h) 


ir(d  +  2A)5 3 

and  the  remaining  term  with  a  little  more  work  yields 


.tvl  —  x2  +  sin  'j 


1 

sa  Id 


(14) 


o(s)  = 


w(d  +  2/2)5" 


c^l  -  x2  (6x2  +  1)- 


sin 


b 

sa  Id 


(15) 


where  b  =  s  Id  for  5  <  d  and  b  =  1  for  5  >  a,  so  the  upper  limit  has  a  branch  point  at  s  =  d 
and  the  lower  has  one  at  s  =  h .  It  is  useful  to  define  the  auxiliary  function 


RA(f,x)  =  Wl  -  x2  f{x)  -  sin  ](jr) 


(16) 


Then  combining  (14)  and  (15),  one  has 


d 

7T (d  +  2/l  )5  3 


d2RA(6x2  + 


U) 


b 

sa  Id 


Ah2H(s 


h)RA(-\,x) 


1 

sa  Id 


The  branch  structure  is  the  same  as  for  c(:  the  lower  limit  contributes  nothing  when  s  <  h .  Furth¬ 
ermore  the  5  <  h  term,  the  upper  limit,  is  part  of  the  total  for  5  >  h .  Both  branches,  s  <  d  and 
5  >  d,  have  this  same  structure.  This  structure  is  utilized  in  the  evaluation,  program  K3CRDN, 
listed  in  the  appendix.  The  programming  is  straight  forward  except  that  this  form  of  the  statement 
function  RX  assumes  SQRT(0.)  returns  zero.  Also  use  of  double  precision  has  been  made  to  avoid 
round-off  error  near  5=0. 

This  routine  checks  normalization  and  average  chord  length  also.  The  probability  (or  integral) 
distribution  is  defined 


C(5)  =  f  c(x)  dx 

J.V 


(18) 


and  the  normalization  condition  is  C(0)  =  1.  Furthermore  there  is  a  theorem,  sometimes  known  as 
the  Cauchy  theorem  [5],  that  three  dimensional  chord  distributions  of  convex  bodies  have  an  average 
chord  length,  5  =4 V /S,  where  V  is  the  volume  and  5,  the  surface  area.  For  a  right  circular 
cylinder  this  is 


5  =  2  hd /(d  +2/2)  (19) 

Five  runs  of  K3CRDN  are  also  exhibited  in  the  appendix  for  different  values  of  the  aspect  ratio,  hid, 
giving  the  sum  over  the  density  to  check  normalization,  the  average  chord  length  calculated  by 

5  =  (  sc(s)ds  (20) 

Jo 


and  that  obtained  from  (19)  for  comparison.  These  runs  were  all  for  2000  bins  (MS)  but  comparing 
series  of  runs  for  fixed  d  and  h  with  number  of  bins  increasing  exhibited  convergence.  The  accuracy 
obtained  of  3  to  4  figures  is  acceptable  for  a  32  bit  machine  except  in  the  case  for  h  Id  =  10,  where 
the  norm  error  is  0.2  percent.  This  is  possibly  improvable  by  resort  to  double  precision  should  the 
need  arise.  However  neither  limit  d  —  0  or  h  —  0  exists  for  the  final  form  so  one  might  expect 
problems  with  either  very  high  or  very  low  aspect  ratio. 


A  slightly  modified  version  of  K3CRDN  has  been  run  to  produce  the  graphs  in  Fig.  1,  showing 
results  for  the  three  most  moderate  aspect  ratios:  0.5,  1.0,  and  2.0.  Note  that  the  jump,  in  both 
cases  where  it  is  isolated  (the  case  h  -  d  =  1  may  be  considered  degenerate  in  that  the  cusp  absorbs 
the  jump),  is  associated  with  the  chords  =  h,  while  the  cusp  is  associated  with  ,v  =  d.  The  values 
of  h  and  d  were  set  at  2/^5  and  4/^5  for  the  high  and  low  aspect  ratios  and  ^ 2  in  the  other  case. 
All  three  cases  have  Max(s)  -  2. 

The  jump  and  cusp  are  interesting  features  from  several  standpoints.  The  differential  chord  dis¬ 
tribution  for  a  rectangular  parallelopiped  (RPP)  also  exhibits  jumps,  three  if  all  three  dimensions  arc 
different  |4|.  Having  observed  these  features  one  can  propound  reasons  for  their  occurrence,  care¬ 
fully  couching  them  in  the  language  of  geometric  probability  |6)  to  mask  any  suggestion  of  hand  wav¬ 
ing.  The  basic  uniform  measure  for  chord  distributions  in  three  dimensions  is  four  dimensional. 
These  four  dimensions  correspond  to  points  on  any  plane  plus  the  element  of  solid  angle  at  each  such 
point.  In  this  qualitative  argument,  one  need  only  point  out  that  when  a  particular  chord  length  has 
an  extra  dimension  available  to  it,  that  length  must  occur  with  higher  probability  in  the  distribution. 
In  the  case  of  the  RPP  one  sees  that  the  occurrence  of  opposite  parallel  faces  provides  this  extra 
dimension  for  that  length  equal  to  the  distance  between  the  two  faces.  Moreover  one  sees  that  this 
length  may  be  approached  continuously  from  above  by  variation  of  solid  angle  but  not  from  below. 
That  is,  a  chord  produced  by  displacing  a  perpendicular  chord  by  a  small  angle  always  makes  a 
longer  chord.  In  other  words  this  feature  indeed  has  the  characteristic  necessary  to  provide  an  “up" 
jump  in  the  chord  distribution.  In  the  case  of  the  circular  cylinder  one  sees  that  “opposing  parallel 
faces”  must  be  expanded  to  the  more  inclusive  term,  “opposing  parallel  elements,"  remembering  the 
elements  of  a  cylinder  are  the  straight  lines  constituting  the  cylindrical  surface.  Furthermore  one  sees 
that  there  is  a  continuous  approach  from  below  as  well  as  from  above.  The  length  of  a  diameter  is 
approached  from  below  by  angular  displacements  perpendicular  to  the  elements  and  from  above  by 
displacements  parallel  to  the  elements.  Hence  opposing  parallel  elements  have  a  cusp  signature  in  the 
chord  distribution. 

The  cusp  and  jump  are  also  of  interest  for,  shall  we  say.  their  mildly  singular  character.  The 
cusp  is  really  singular  but  “mildly"  so,  as  it  is  integrable.  Having  seen  above  that  the  cause  of  these 
features  is  opposing  parallel  elements,  one  also  sees  that  such  features  are  pervasive  among  typical 
sensitive  volumes  in  memory  cells.  While  a  depletion  region  (that  is,  a  sensitive  volume)  is  defined 
in  part  electrostatically  and  thus  must  have  rounded  corners,  such  rounding  is  likely  to  be  mainly 
local,  leaving  still  significant  areas  of  opposing  parallel  elements.  Indeed  while  this  electrostatic- 
effect  tends  to  ameliorate  parallel  faces,  it  may  enhance  cylindrical  parallelism.  It  seems  there  is  rea¬ 
son  to  expect  physical  consequences  from  these  features.  In  the  evaluation  of  the  upset  rate  in  space 
the  coincidence  of  one  of  these  mildly  singular  features  with  a  sudden  jump  in  cosmic  ray  flucnce  at  a 
particular  LET  should  provide  an  especially  dramatic  onset.  Of  itself  this  feature  is  not  qualitatively 
distinguishable  because  already  jumps  in  the  LET  spectrum  occur.  However  one  prominent  evalua¬ 
tion  method,  Monte  Carlo,  would  miss  enhanced  onset  from  this  source  because  it  is  incapable  of 
evaluating  mildly  singular  features. 

Some  limit  considerations  influenced  the  evaluation  while  others  did  not.  It  seems  worthwhile  to 
mention  them  for  their  possible  affect  on  modifications.  The  jump  is  finite  and  is  no  problem.  No 
provision  was  made  for  the  case  that  a  bin  center  point  might  fall  on  the  cusp,  which  is  a  singularity. 
This  point  is  simply  not  evaluated  in  the  direct  coding.  However  in  the  evaluation  of  the  two  elliptic- 
integral  library  functions,  a  stop  occurs  if  both  the  first  two  arguments  are  zero.  This  could  only  hap¬ 
pen  if  a  bin  center  fell  on  a  cusp  within  computer  precision.  A  weak  attempt  to  avoid  this  was  made 
by  dividing  bins  irrationally  with  respect  to  s  =  d  or  h  by  using  maximum  chord  length.  Since  this 
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would  have  imposed  unusual  input  restrictions,  cases  of  rational  ratios  d  or  h  to  Max(s)  have  been 
evaluated  and  no  problem  was  encountered.  The  case  where  one  formed  bins  by  dividing  d  or  h 
instead  of  Max(s)  has  not  been  investigated.  In  fact  an  investigation  of  the  .v  —  d  limit  when  h  =  d 
was  made  using  known  limit  properties  of  the  Legendre  forms  of  elliptic  integrals.  As  indicated  ear¬ 
lier  this  contributed  little  to  the  results  except,  perhaps,  a  familiarity  with  problems  avoided.  The  dis¬ 
tribution  end  points  are  easily  evaluated.  If  s~  =  dr  +  /t2, 

c(s„,)  =  0  (21) 

c(0)  =  16/(3ir(d  +  2h ))  (22) 

the  latter  arising  as  a  limit  of  the  c0  term.  Thus  evaluation  near  s  =  0  should  be  avoided  without 
making  special  provision.  In  Fig.  1,  the  three  graphs,  in  order  of  increasing  aspect  ratio,  begin  at 
0.474,  0.400,  and  0.380,  respectively,  according  to  (22). 

Perhaps  it  should  be  remarked  that  the  rigorously  proper  mathematical  milieu  for  some  of  the 
manipulations  above  is  non-standard  analysis  since  for  completeness  delta  functions  are  needed.  That 
it  was  nearly  possible  to  carry  through  this  analysis  without  reference  to  delta  functions  is  interesting 
in  itself.  The  one  reference,  page  2  and  1,  was  only  to  point  out  that  by  using  one,  a  difficulty  can 
be  seen  not  to  be  a  problem.  Probably  the  reason  this  short  cut  was  feasible  is  that  the  delta  function 
fortuitously  integrates  to  zero.  This  is  proven  by  the  following  exercise. 

The  same  methods  as  above  were  applied  to  evaluate  Kellerer’s  integral  chord  distribution  for  a 
circular  cylinder.  An  ambiguity  arose  as  to  the  constant  of  integration  applicable  for  the  s  <  d 
branch  which  is  related  to  the  choice  of  the  use  of  arccosine  or  arcsine  in  the  evaluation  of  one  of  the 
standard  integrals.  This  ambiguity  was  easily  resolved  by  enforcing  the  normalization  condition,  fol¬ 
lowing  equation  18.  At  the  same  time  this  results  in  continuity  at  s  =  d .  This  proves  the  delta  func¬ 
tion  makes  no  contribution  because  of  contributing  one  would  show  as  a  jump  at  s  =  d .  No  ambi 
guity  arises  for  the  ,v  >  d  branch  because  it  approaches  zero  at  maximum  s  due  to  the  coalescence  of 
the  limits  of  integration. 

The  integral  distribution,  evaluated  directly,  agreed  entirely  with  the  numerical  integral  of  the 
differential  distribution  thus  providing  an  additional  confirmation  of  these  evaluations.  The  direct 
integral  is  available  but  is  not  included  here  because  integrating  the  differential  is  just  as  efficient 
whenever  the  complete  distribution  is  being  evaluated.  Furthermore  the  integral  algorithm  seems  to 
be  slightly  more  sensitive  to  the  need  for  higher  precision  at  small  x . 

This  same  evaluation  using  Legendre  elliptic  integrals  has  been  published  by  U.  Madcr  |7) 
along  with  a  number  of  related  results.  For  example,  the  infinite  limit  of  both  d  and  h  are  obtained 
and  the  latter  gives  our  c ,  term.  Graphical  results  appear  to  be  the  same  after  renormalizing  to 
account  for  the  use  of  different  dimensions.  The  author  wishes  to  thank  A  M.  Kellerer  for  this  refer- 
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$  R  K3CRDN 

ENTER  D , H , MS , KS , FO : 

2. ,2. , 2000 , 0 , 0 

NORM:  0.9999697  AVE :  1.333315  COMPARE:  1.333333 

$  R  K3CRDN 

ENTER  D , H , MS , KS , FO : 

2 .  , 1 .  ,2000 ,0,0 

NORM:  1.0005059  AVE:  0.999779  COMPARE:  1.000000 

$  R  K3CRDN 

ENTER  D , H , MS , KS , FO : 

1 .  , 2 .  , 2000 , 0 , 0 

NORM:  1.0005429  AVE:  0.800043  COMPARE:  0.800000 

$  R  K3CRDN 

ENTER  D , H , MS , KS , FO : 

2 . , . 2 , 2000 , 0 , 0 

NORM:  1.0000465  AVE:  0.333263  COMPARE:  0.333333 

$  R  K3CRDN 

ENTER  D , H , MS , KS , FO : 

.2,2. ,2000,0,0 

NORM:  0.9978931  AVE:  0.190463  COMPARE:  0.190476 


